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ABSTRACT 

The z = 0.9902 system along the quasar PG 1634 + 706 line of sight is a strong 
Mgll absorber (1^^(2796) > 0.3 A) with only weak Civ absorption (it is "Civ- 
deficient"). To study this system, we used high-resolution spectra from both the 
Hubble Space Telescope (HST)/Space Telescope Imaging Spectrograph (STIS) 
and the Keck I telescope/High Resolution Echelle Spectrometer (HIRES). The 
STIS spectrum has a resolution of i? = 30, 000 and covers key transitions, such 
as Sill, Cll, Sim, Cm, Siiv, and Civ. The HIRES spectrum, with a resolution 
of i? = 45,000, covers the Mgl, Mgll and Fell transitions. Assuming a Haardt 
and Madau extragalactic background spectrum, we modeled the system with a 
combination of photoionization and collisional ionization. Based on a comparison 
of synthetic spectra to the data profiles, we infer the existence of the following 
four phases of gas: 

i) Seven Mgll clouds have sizes of 1-1000 pc and densities of 0.002-0.1 cm~^, 
with a gradual decrease in density from blue to red. The Mgll phase gives rise to 
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most of the C IV absorption and resembles the warm, ionized inter-cloud medium 
of the Milky Way; 

ii) Instead of arising in the same phase as Mgll, Mgl is produced in separate, 
narrow components with b ~ 0.75 km s~^. These small Mgl pockets (~ 100 AU) 
could represent a denser phase (~ 200 cm~^) of the interstellar medium (ISM), 
analogous to the small-scale structure observed in the Milky Way ISM; 

iii) A "broad phase" with a Dopplcr parameter, 6 ~ 60 km s~^, is required 
to consistently fit Lya, Ly/3, and the higher-order Lyman-series lines. A low 
metaUicity (logZ < —2) for this phase could explain why the system is "ClV- 
deficient", and also why Nv and Ovi are not detected. This phase may be a 
galactic halo or it could represent a diffuse medium in an early-type galaxy; 

iv) The strong absorption in Si IV relative to CiV could be produced in an 
extra, collisionally ionized phase with a temperature of T ~ 60, 000 K. The col- 
lisional phase could exist in cooling layers that are shock-heated by supernovae- 
related processes. 

Subject headings: quasars — absorption lines; galaxies — evolution; galaxies — 
halos 



1. Introduction 

Quasar absorption lines provide a unique and intriguing way to determine the kine- 
matics, chemical composition, and ionization state of gas in intervening galaxies. For most 
galaxies, a multi-phase interstellar medium (ISM) (i.e. a medium with different densities 
and temperatures) is likely to exist (McKee & Ostriker 1977; Giroux, Sutherland, & ShuU 
1994). Thus, it follows that any random line of sight passing through a galaxy is hkely to 
encounter several phases of gas. 

Strong Mgll absorbers (Wr(2796) > 0.3 A) are almost always found within an impact 
parameter of 35h~^ kpc of a luminous galaxy (> 0.05L*), where L* is the Schechter luminosity 
(Bergeron & Boisse 1991; Bergeron et al. 1992; Le Brun et al. 1993; Stcidcl, Dickinson, & 
Persson 1994; Steidel 1995; Steidcl et al. 1997). In addition, the strong CiV absorption that 
is characteristic of most such absorbers is thought (Churchill et al. 1999, 2000) to arise in 
a corona similar to that surrounding the Milky Way disk (Savage et al. 1997). 

We study the z = 0.9902 system along the line of sight to the quasar PG 1634 -|- 706 
{zem = 1.335). This system is a strong Mgll absorber with at least five blended components 
and with a deficiency of CiV relative to other systems with similar Wr{Mgu) (Churchill 
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et al. 2000). In a previous study, Charlton et al. (2000) derived constraints on the physical 
properties of the system, based upon the combination of low-resolution data from the Faint 
Object Spectrograph (FOS) on-board the Hubble Space Telescope (HST) and high-resolution 
data from the High Resolution Echcllc Spectrometer (HIRES) on the Keck 1 telescope. At 
that time, only low-resolution data were available for many key transitions, including low- 
ionization tracers (Sill and Cll), high-ionization tracers (Si IV and Civ), and the Lyman- 
series hues. However, the HIRES profiles of Mgl, Mgll, and Fell were used to form a template 
that was applied to model the low-resolution data. 

Based on the FOS and HIRES data, Charlton ct al. (2000) came to the following 
conclusions: (1) The mctallicity (expressed in the solar units) of the Mgll phase would 
be low (—1.5 < logZ < —1) if it was to produce a full Lyman hmit break; (2) A broad 
component (40 < 6(H) < 80 km s~^) was needed to fit the Lya profile; (3) An additional 
broad component also appeared to be needed to produce CiV; (4) Si IV could not fully arise 
in the inferred low- or high-ionization phases, which suggested the existence of a blend or 
an additional component. 

However, some important issues still remained unresolved: (1) whether the Mgl could 
be produced in the same clouds as Mgll; (2) whether the ionization conditions vary across 
the profiles; (3) whether CiV is produced by a single, broad phase or by several separate, 
narrower clouds; and (4) whether Si IV is affected by a blend and, if not, how is it produced? 

The first issue is important because Mgl is under-produced by the Mgll phase in many 
strong Mgll absorption systems (Churchill 1997; Churchill, Vogt, & Charlton 2002). In ad- 
dition, Mgl absorption is often very strong in the environment where damped Lya absorbers 
(DLAs) arise and therefore could be used to trace the origin of the strongest Hi absorption. 
A resolution to the Mgl/Mgll issue in the case of the z = 0.9902 Mgll system could serve 
as a basis for solving these more general problems. 

The third issue is of particular interest because this system is Civ-deficient (Churchill 
et al. 2000). Most strong Mgll absorbers have comparable strong CiV absorption, which 
is likely to arise in a corona similar to that surrounding the Milky Way. Although CiV is 
detected in the low-resolution FOS spectrum, it is relatively weak. Therefore, the additional 
component suggested by Charlton et al. (2000) may not be analogous to a corona. The 
CiV profiles could be resolved to determine whether the "corona phase" is absent, or just 
weak in CiV. Either low metallicity or high-ionization conditions could lead to weak CiV 
absorption from a corona. 

With the release of grating E230M spectra from the Space Telescope Imaging Spectro- 
graph (STIS) on-board the HST, we now have high-resolution coverage of the Lyman series. 
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Sill, Cll, Sim, Cm, SilV, and Civ transitions. In this paper, we present the results from 
photo-/coUisional ionization modeling of the z ~ 0.9902 Mgll system. Combining the newer 
data with those on the Mgl, Mgll, and Fell from Keck/HIRES, we determine the minimum 
number of phases required to produce the observed absorption in the many detected chemi- 
cal transitions. For each phase, we place constraints on physical properties such as densities, 
temperatures, and sizes of the various environments that give rise to the absorption features 
displayed in the observed spectral profiles. 

In § 2, we briefly describe the data that we use to constrain our models. A summary 
of OTU modeling techniques is presented in § 3 and our major results are outlined in § 4. 
Finally, in § 5, we summarize and give a physical interpretation of the system. 

2. Data and Analysis 

2.1. Keck/HIRES 

The optical part of the spectrum (with Mgl 2853, Mgll 2796, Mgll 2803 and Fell 2600 
covered between 3723 A and 6186 A) was obtained from the HIRES on the Keck I telescope 
on July 4 and 5, 1995 (Vogt et al. 1995). This spectrum has a resolution oi R — 45,000 
(FWHM ~ 6.6 km s""^) and a signal-to-noise ratio of ~ 50 per three-pixel resolution element 
(Charlton et al. 2000). 

The HIRES spectrum was reduced with the IRAF APEXTRACT package for echelle 
data and was extracted using the optimal extraction routine of Horne (1986) and Marsh 
(1989). The wavelengths were cahbrated to vacuum using the IRAF task ECIDENTIFY, 
and shifted to heliocentric velocities. Continuum fitting was based on the formalism of 
Sembach & Savage (1992). 

Voigt profiles were used to simultaneously fit Mgl, Mgll, and Fell, with free parameters 
redshift, column density and Doppler parameter, 6, for each component. The program 
MINFIT (Churchill 1997), using a formalism, found the minimum number of components 
required to fit the system. The Mgl 2853, Mgll AA2796, 2803 and Fell 2600 transitions, 
covered by HIRES, are shown in velocity space in Figure 1. Other Fell transitions were also 
covered, but Fell 2600 has the highest signal-to-noise. 

The fitting procedure was summarized in Churchill & Vogt (2001), and described in 
detail in Churchill (1997). Components are dropped from the fit unless retaining them 
produces a value that is significantly better (97 % by the F-test). Errors in redshifts, 
column densities, and Doppler parameters were determined, using the diagonal elements of 
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the covariance matrix, determined from the curvature matrix. Components were rejected if 
their total fractional errors exceeded 1.5. 



2.2. HST/FOS 

Before the release of the STIS data, low-resolution UV spectra were obtained using 
the G190H and G270H gratings of HST/FOS. These spectra were presented in Impey et al. 
(1996) and Bahcall et al. (1996). Since the higher-resolution STIS spectrum covers most of 
the key transitions for the z = 0.9902 system, we only used the FOS spectra to study 
the Lyman limit break at 1830 A, which implies the neutral hydrogen column density 
logA^(Hl) > 18. The data reduction, wavelength calibration, fluxing, and continuum fit- 
ting of the FOS spectra were performed as part of the QSO Absorption Line Key Project 
(Bahcall et al. 1996). 

2.3. HST/STIS 

Two high-resolution (/? = 30, 000) datascts, with different wavelength coverage, were 
obtained with HST/STIS and retrieved from the HST data archive. The E230M spectrum, 
obtained by S. Buries (May and June 2000), has a wavelength coverage of 1865 A to 2673 A, 
and a total exposure time of 29, 000 seconds. The E230M spectrum, provided by B. T. 
Jannuzi (June 2000), has a wavelength coverage of 2303 A to 3111 A, and a total exposure 
time of 26,435 seconds. The Lyman series. Sill, Cll, Sim, Cm, SilV, and CiV are all 
covered by the E230M STIS spectra. A G230M spectrum, covering 1830 A to 1870 A, was 
also obtained by S. Buries (May 2000), but it did not cover any key transitions for our 
analysis of the system. 

The STIS spectra were reduced using the standard pipeline (Brown et al. 2002). The 
extracted spectra were averaged between different exposures and the continuum fitting was 
performed with standard IRAF tasks (Churchill & Vogt 2001). 

The spectra of various transitions are shown in Figure 1. Ly7 is apparently blended 
with another feature, based on its asymmetry relative to the other Lyman-series transitions. 
LyS is not shown due to its noisy spectrum and its contamination by a blend. Sill 1190 is 
blended with Nv 1243 from the z = 0.9056 system. Cll 1335 is contaminated by SilV 1394 
from the z = 0.9056 system, but based on SilV 1403, the contribution from SilV 1394 would 
be negligible. A relatively strong, unidentified feature is apparent to the red of the Cll 1335 
profile. Sim 1207 is blended with Sill 1260 from the z = 0.9056 system. Only Nv 1239 
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is shown because no absorption is detected in either of the Nv AA1238, 1242 profiles and 
Nv 1239 is the stronger of the two. Ovi 1032 is not shown. It is clearly inconsistent with 
the undetected Ovi 1038. 



3. Methods for Modeling 

The goal of the present study was to determine the physical conditions under which the 
spectra shown in Figure 1 were produced. To begin, we applied a Voigt profile fit to the 
Mgll transitions to determine the locations, column densities, and Doppler parameters of 
the components required to reproduce the observed absorption. We assumed that each one 
of these components was produced by an individual "cloud" and we modeled each one of 
these clouds with the photoionization code Cloudy, version 90.4 (Ferland 1996). 

The clouds were assumed to be constant-density, plane-parallel slabs in photoionization 
equilibrium. They were ionized by a 2; = 1 Haardt and Madau (1996) extragalactic back- 
ground spectrum normalizing with an ionizing photon density of logn^ = —5.2. Alternative 
input spectra were also briefly explored, as discussed in § 4.5. To simulate these clouds, the 
required input parameters were: i) the column density of Mgll, iV(Mgll), ii) the metallicity, 
Z (expressed in the units of the solar value), iii) the ionization parameter, U (defined as 
the ratio of the number density of photons capable of ionizing hydrogen, n^^ to the number 
density of hydrogen, uh)-, and iv) the abundance pattern. The column density of Mgll was 
obtained from the Voigt profile fit. Both the metallicity and the ionization parameter were 
chosen from a range reasonable for interstellar and intergalactic clouds (i.e. —3 < log ^ ^ 
and —5 < logf/ < 0), with the initial assumption of a solar abundance pattern. If a si- 
multaneous fit to all transitions was not possible with solar abundances, alternatives were 
explored. 

Once the input parameters were specified, we ran Cloudy to optimize on the column 
density of a selected species, such as Mgll. During the run. Cloudy used an initial guess of 
iV(Hl) and calculated the corresponding A^(Mgll). The output A^(Mgll) was compared to the 
value measured from the Voigt profile fit, and the value of A^(Hl) was adjusted accordingly 
for the next iteration. This process was repeated until the difference between the output 
A^(Mgll) and the designated one was neghgible. 

Upon the completion of the final iteration, the following parameters were retrieved for 
the resulting cloud: the column density of every species with spectral coverage, the size 
of the cloud, and the temperature. The temperature was used to calculate the thermal 
component of the Doppler parameter of each ion, according to fct^jgrm ~ 2kT/m, where m 
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is the atomic weight of the ion. The turbulent component of b, which is due to the bulk 
motion of the gas and is the same for all the ions, was calculated from the measured 6(Mg) 
using 6^^^^ = 6^(Mg) — b^^ermi^s)- Combining the turbulent and thermal components, we 
obtained the Doppler parameter of all other individual ions. 

Cloudy was run for all five Mgll clouds, each time calculating and b for each species. 
A "pure" spectrum was then produced for each transition from the N and b of all five 
components. Then, we generated a synthetic spectrum by convolving the "pure" spectrum 
with the instrumental spread function. Finally, the synthetic spectrum was superimposed 
on the observed one for comparison. 

We experimented with various metallicities and ionization parameters and attempted to 
obtain a synthetic spectrum that displayed minimal discrepancy with the observed one. A 
visual comparison of goodness of fit was used to reject models, as illustrated in Charlton et al. 
(2000) and Rigby et al. (2002). This was adequate to confidently constrain parameters to a 
relatively narrow range within the very large parameter space. The log U and log Z values 
given below are accurate to about 0.1 dex. We considered using a quantitative approach, 
but found it to be impractical because a single pixel or group of pixels for a particular 
transition often dominated the assessment. 

We encountered several issues that were unresolvable using a single phase of gas: Mgl 
was under-produced for any choice of U; the shape of the Lya profile could not be fit; the 
velocity centroids of Si IV and CiV were offset from that of Mgll; and the Si IV and CiV 
transitions could not be fit simultaneously. Additional clouds were included in the Mgll 
phase to account for the offset Siiv and Civ absorption. The resolution of the other three 
issues requires the inclusion of independent phases: A colder, denser phase is needed to 
produce the observed absorption in Mgl. A broad, low-metallicity component is needed to 
account for the shape of the Lja profile. Finally, a collisionally ionized phase is needed to 
fit Si IV without over-producing CiV. Each of these phases will be explained in detail in § 4. 

4. Results 

Here we present constraints on each of the four phases required to fit the z = 0.9902 
system, beginning with the lowest-ionization phase. Table 1 provides a summary of the 
cloud properties for all four phases in a plausible model that satisfies all constraints. 
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4.1. Mgl Phcise 

The five Mgll clouds from the original Voigt profile fit could not produce the ob- 
served equivalent width of Mgl, regardless of the choice of ionization parameter. The ratio 
A^(Mgl)/iV(Mgll) is nearly constant at ~ 0.01 over a large range of ionization parameters, 
U, and for the range of metallicity appropriate for this system. 

In order to obtain the observed equivalent width ratio Wr {2853} /Wr{2796) , a difference 
in the curve of growth between the two transitions was used, as shown in Figure 2. In general, 
when an absorption profile is unsaturated, its equivalent width, Wr, increases linearly with 
its column density, N (i.e. on the linear part of curve of growth). However, when the 
absorption becomes saturated, Wr stays almost constant with N, at a value that depends 
only on Doppler parameter, b (i.e. on the flat part of curve of growth). In Figure 2, our 
typical iV(Mgll) is on the fiat part of the curve of growth, while A'"(Mgl) is on the linear 
part. Thus, we could produce a larger Wr{28b3)/Wr{2796) using a smaller 6, in a "Mgl 
phase" . 

A unique Voigt profile fit to Mgl could not be obtained. We assumed five clouds 
at the same velocities as the broader Mgll cloud components, with a constant Doppler 
parameter and a constant ionization parameter. The Doppler parameter is constrained to 
be 6 < 1 km s^^, in order that Mgll is not over-produced by the combination of this Mgl 
phase and the Mgll cloud phase (see § 4.2). However, the ionization parameter has to 
be logC/ < —8 for 6 = 1 km s~^, in order that Ol is not over-produced. Alternatively, 
^ogU < —7.5 if 6 = 0.75 km s~^ (Ol is on the fiat part of the curve of growth). Column 
densities for a sample model with 6(Mg) = 0.75 km s~^ and log U — —7.5 are listed in Table 
1. With 6(Mg) = 0.75 km s"^, the temperature is constrained to be lower than 800 K. Cloudy 
yields an effective temperature close to this value for the clouds in our sample model. An 
ionization parameter of logf/ = —7.5 corresponds to a cloud density of n{H) ~ 200 cm~^. 
The cloud sizes are quite small, ranging from 60-160 AU. 

If the b parameter is small, the Mgl phase will have trivial contribution to the overall 

equivalent width of any high-ionization metal transition. In Figure 3, the dotted lines repre- 
sent the contribution from the Mgl phase for the model listed in Table 1. The contribution 
to Mgll is negligible, therefore Mgll has to be entirely produced by the Mgll cloud phase 
(see § 4.2). 

Assuming the five Mgl clouds have the same metallicity, the metallicity of this phase 
is constrained to be logZ > —0.8, by the Lya profile. At logZ = —0.8, the Hi column 
densities of the five Mgl clouds are in the range of 17.3 < logA^(Hl) < 17.7, as listed in 
Table 1, yielding an effective total of logiV(Hl) ~ 18.1. Thus, at this metallicity, the Mgl 
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phase will give rise to the full Lyman limit break detected in the FOS spectrum, which 
requires logA'"(Hl) > 18. 

4.2. Mgii Phase 

The column densities and b parameters of the five Mgll clouds were obtained from a 
Voigt profile fit to the residuals in Mgll unaccounted for by the Mgl phase. The five clouds 
in this phase have considerably larger b parameters (6 < b{Mg) < 15 km s~^) than the Mgl 
clouds, as listed in Table 1. 

With the detection of strong Fell in clouds Mgll^^, Mgllg, and Mgllg, the ionization 
parameters of these clouds are constrained by A^(Fell)/A^(Mgll). The remaining two clouds 
(Mgll4 and Mgllg) must have higher ionization parameters, since Fell is very weak. Si IV was 
used to determine their ionization states, with the assumption that the Si IV transitions arise 
in the same phase as Mgll. Thus, the five clouds have ionization parameters of log U ~ —3.7, 
—4.0, —3.5, —3.0, and —2.5, respectively, as listed in Table 1. 

The Si IV and CiV profiles arc not aligned with the bulk of the absorption seen in Mgl, 
Mgll, and the other low-ionization transitions. Thus, additional offset clouds with higher 
ionization parameters are needed. A Voigt profile fit to Siiv yielded two additional clouds 
(Mgllg and Mgll^) at f ~ 42 and f ~ 55 km s^"*^, respectively. Again, column densities and 
b parameters of these two clouds are listed in Table 1 . Ionization parameters are constrained 
to be logU ~ -2.5 and -2.7, by A^(Silv)/A^(Clv). This is also consistent with the weak 
absorption displayed in the red wings of the Mgll and Sill profiles. 

A wide range of ionization states (—4.0 < logf/ < —2.5) is displayed in this Mgll phase, 
with a tendency for increase in ionization stage from blue to red in the spectrum. Therefore, 
most of the absorption in both low-ionization transitions, such as Mgll, Fell, Sill, and Cll, 
and high-ionization transitions, such as Sim, Cm, and CiV, is produced in this phase, as 
shown in Figure 3. However, the substantial under-production of Si IV suggests the existence 
of an extra component, as discussed in § 4.4. 

Assuming that all seven clouds in the Mgll phase have the same metallicity, the metal- 
licity is constrained to be logZ > —0.65 by the Lyman series. However, even at the lower 
hmit of logZ = -0.65, with 14.9 < log A'"(Hl) < 17.1, the Mgll phase contributes insignif- 
icantly to the observed full Lyman limit break. Alternatively, if the Mgl and Mgll phases 
are assumed to have the same metallicities, a metallicity of \ogZ ~ —0.8 is required by the 
Lya profile (Figure 4). As mentioned in § 4.1, at this metallicity the five Mgl clouds will 
give rise to the observed full Lyman limit break (see Table 1). 
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4.3. Broad Hi Phase 

The Lja profile is too broad to be fully produced by the Mgl and Mgll clouds, as 
is shown in Figure 4. However, if the metallicity was decreased to give rise to more Lya, 
the narrower Lyy would be over-produced. The curve of growth of the Lyman-series lines 
(Churchill & Charlton 1999) allows for a centered, broad component that can resolve this 
dilemma (i.e. it accounts for the broader, deeper Lya profile, as well as the wing structure 
in Ly/3, without over-producing the narrower Ly7 profile.) 

A Voigt profile fit to Lya yielded the column density and b parameter of the broad phase. 
If 6(H) > 70 km s~^ and A^(Hl) < 10^^ cm~^, we can produce the saturated absorption 
from —100 to +120 km s~^, but the wings will be over-produced; if 6(H) < 50 km s~^ 
and A'^(Hl) > 10^^ cm~^, we can produce the range of the saturated absorption, but the 
wings will be under-produced. The best fit was found to be with A^(Hl) ~ 10^^'^ cm~^ and 
6(H) - 58 km s'K 

In principle, the broad component could be either photoionized or coUisionally ionized. 
In both cases, the metallicity is constrained so as not to over-produce any of the high- 
ionization transitions (i.e. CiV, Nv and Ovi). The Nv and Ovi are not detected, and 
Civ is fully produced in the Mgll phase (as required by the structure appearing in the CiV 
profiles.). 

For photoionization, an upper limit of log Z < —2 is placed on the metallicity. regardless 
of the choice of ionization parameter. Otherwise, either low- or high-ionization transitions 
would be over-produced. If the metallicity is \ogZ < —2.5, there will be no constraint on 
the ionization parameter from the metal transitions. None of them is significantly produced 
at such a low metallicity. However, a rough upper limit of logU < — 1 can be estimated by 
the requirement of a reasonable cloud size (less than several hundred kpcs). A good fit for 
this phase was obtained for a model with a metallicity of log Z ~ —2.5 and an ionization 
parameter of log U ~ —1.5, as fisted in Table 1. This model gives a cloud ~ 60 kpc in size. 
The fit to the Lyman series is shown by the dashed-dotted curve in Figure 4. 

For coUisional ionization with pure thermal broadening (i.e. no bulk motion), the tem- 
perature would be T ~ 10^-^ K, using 6^ = 2kT/m for Hi. At this temperature, the 
metallicity has to be exceptionally low, logZ < —3.5, in order not to over-produce the 
high-ionization transitions. If, instead, bulk motion contributes, the b parameter will be 
larger for each metal species than the one given by pure thermal scaling of 6(H). For ex- 
ample, as the temperature decreases from T ~ 10^ '^ K, A^(Clv) increases towards a peak at 
T ~ 10^ *^ K. This, together with a larger 6(C), gives a larger Civ equivalent width when bulk 
motion is involved. In order not to over-produce CiV, the metallicity has to be even lower 
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than the pure thermal case. A metaUicity as low as logZ ~ —3.5 is unlikely. Therefore, 
coUisional ionization is ruled out for this phase. The broad Hi absorption is more likely to 
arise in a photoionized cloud, but even then the metaUicity is constrained to be quite low 
(logZ<-2). 

4.4. CoUisional Ionization Phase 

The broad, smooth absorption profiles in the Si IV transitions are not consistent with 
the structure apparent in the Civ profiles. Because Siiv is offset relative to Lja, it also 
cannot be produced by the centered, broad Hi phase. The need to produce Si IV without 
contributing a substantial amount of absorption to other transitions suggests the existence of 
a coUisionally ionized phase. A Voigt profile fit to the residuals in Si IV yielded a component at 
V ~ 17 km s^^, with A^(Silv) ~ 10^^'^ cm"^ and 6(Si) ~ 17 km s^^. However, for coUisional 
ionization. Si IV peaks at logT ~ 4.8 (Sutherland & Dopita 1993), which corresponds to 
hherm ~ 6 km s~^. Therefore, bulk motion must also contribute to the overall velocity 
dispersion. The temperature range was found to be 10^-^ T < 10^-^ K, in order not to 
over-produce the transitions of higher- and lower-ionization stages. The metaUicity of this 
phase could be set to the same value as the Mgl and Mgll phases, logZ ~ —0.8, but it is 
not well constrained. A lower limit of log Z > —3 is placed by the the Lyman-series lines. 

4.5. Effects of Assumed Input Spectrum 

For all the previous results, we have used the z — 1 Haardt and Madau extragalactic 
background radiation (EBR) spectrum (Haardt & Madau 1996) as the only source of ionizing 
photons. Here we consider the likelihood of additional stellar contributions and consider the 
effect of such alternative spectra in a couple of extreme examples. 

At redshift z ~ 1 the extragalactic background radiation is more likely to dominate 
stellar contributions than it is at other redshifts. At lower redshifts, the amplitude of the 
Haardt and Madau extragalactic background is lower (by a factor of seven at z = (Haardt 
& Madau 1996)). At much higher redshifts, starbursts are common and the amplitude of 
the EBR from quasars has levelled off. 

Unfortunately, in the case of the PC 1634 + 459 line of sight, we have no information 
about the colors or morphologies of the host galaxies. Therefore, we cannot rule out a 
starburst host for the z = 0.9902 absorber. Only an extreme starburst would affect the 
results, and only within ~ 6 kpc of the center of such a galaxy would the stellar contribution 
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strongly dominate the EBR (1% of the photon flux of 10^^ would escape the burst region 
(Hurwitz, Jelinsky, & Dixon 1997)). We would expect the absorption from the central 
region of a starburst to be stronger than it is for this z = 0.9902 system (Bond et al. 2001). 
Furthermore, we would expect to pass through some layers of highly ionized gas, and this 
absorber is weak in CiV and undetected in Nv and Ovi. Therefore, a strong starburst host 
galaxy is unlikely, but we still consider here the effect of a couple of modified spectra. 

Many different spectral shapes would be possible once we allow a stellar contribution 
and it is unrealistic to consider them all, so we choose two extreme cases: 1) a 0.01 Gyr 
instantaneous burst, which is characterized by the edges of Hi, Hcl, and Hell; and 2) a 
0.1 Gyr instantaneous burst with an extreme Hi edge, but with a relatively flat spectrum 
above that edge, similar to the Haardt and Madau shape. Both stellar spectra, assuming 
solar metallicity and a Salpeter IMF, were taken from Bruzual & Chariot (1993). We have 
run models with the burst contribution set to ten times that of the EBR, at 1 Rydberg. We 
begin with the model in Table 1 and consider how the addition of the burst would affect our 
conclusions. In a more realistic model, however, the spectral shape would be a function of 
position, with the EBR likely to be dominant at some locations, but the solution to that case 
would be in between the extreme starburst models and a pure Haardt and Madau model. 

Stellar contributions from a non-burst model arc more likely to dominate the EBR at 
energies less than 1 Rydberg. This could affect the ionization balance between the neutral 
and singly ionized transitions. Because of its extreme Lyman limit break (almost a factor of 
1000), the 0.1 Gyr instantaneous burst model can also be used to represent the rough effects 
in this case. 

For the 0.01 Gyr model, due to the softer spectrum (because of the edges of helium), 
the Civ is under-produced compared to the pure Haardt and Madau case. For the same 
reason, for the clouds optimized on Si IV, the lowcr-ionization transitions are found to be 
over-produced compared to Si IV. With this spectral shape, we cannot produce even the 
weak, observed CiV. In principle, this could relate to why this system is CiV-deficient, but 
in practice it is difficult to understand how such a strong stellar flux could affect a large 
enough region to reduce the high-ionization absorption at all velocities. The metallicity 
that we would infer for the 0.01 Gyr model would be quite similar (within 0.1 dex) to that 
for the pure Haardt and Madau case. 

For the 0.1 Gyr model, there are relatively fewer ionizing photons just above the Hi 
edge, and so more hydrogen remains neutral. This gives rise to a stronger Lya line. To 

match the observed Lya profllc. the metallicity would have to be increased, but only by 
~ 0.3 dex. Due to the similar shape to the Haardt and Madau spectrum above the edge, the 
ionization conditions are not very different for metal-line transitions. Slightly more CiV is 
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produced in the 0.1 Gyr model, but this would require an ionization parameter adjustment 
of only ~ 0.1 dex. It is also important to note that for neither burst model is our conclusion 
altered about the need for an additional Mgl phase. However, the 0.1 Gyr burst model 
does give rise to a larger amount of O I in the Mgl clouds, which would somewhat change 
constraints on the cloud densities and sizes. 



5. Summciry and Discussion 

Based upon modeling of high-resolution HST/STIS and Keck/HIRES absorption pro- 
files of the z = 0.9902 strong Mgll absorber toward PG 1634 + 706, we derived the physical 
conditions of the phases of gas encountered along the line of sight. Here, we summarize 
the physical conditions (densities, temperatures, size of structures) that characterize each 
of the phases. These properties indicate relationships between the gas phases and various 
structures in the Milky Way, as well as in other galaxies in the local universe. 

5.1. Mgii Clouds 

Seven blended components with b parameters ranging from 4-15 km s^^ provide a con- 
sistent fit to the Mgll AA2796, 2803 profiles, which extend over ~ 75 km s~^ in velocity. 
Under the assumption of a Haardt and Madau ionizing spectrum, these photoionized clouds 
have densities of 0.002-0.1 cm~^ (corresponding to —4.0 < logU < —2.5), with a gradual 
decrease from blue to red, and temperatures of ~ 10,000 K. CiV was also fully produced 
in this phase. We infer a metaUicity of logZ ~ —0.8, assuming that these clouds have the 
same metallicity as those of a separate Mgl phase. The cloud sizes/thicknesses range from 
parsecs to hundreds of parsecs, as hsted in Table 1. 

The physical conditions in this phase are similar to those in the warm, ionized inter- 
cloud medium of the Milky Way (McKee &: Ostriker 1977). They are also similar to most 
other Mgll clouds in strong Mgll absorbers (Churchill et al. 2000). The similarity with the 
Milky Way disk ISM and the connection between strong Mgll absorbers and the majority 
of L* galaxies suggests that the Mgll absorption comes from the warm ISM of a spiral disk. 
However, strong Mgll absorption can also come from early-type galaxies (Churchill, Steidel, 
& Vogt 1996). In fact, Churchill et al. (2000) found that the three CiV-deficient Mgll 
absorbers in their sample had among the reddest B — K colors and the highest luminosities. 
The similarity of the low-ionization gas properties between "classic" strong Mgll absorbers 
and CiV-deficient Mgll absorbers would then suggest that early-type galaxies house regions 
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with an ISM similar to that in spirals. 

5.2. Pockets of Mg I 

Perhaps the most interesting and surprising result of our modeling was our inference of 
the existence of a cold phase {T < 1000 K). It consists of dense (n(H) ~ 200 cm~^) pockets 
that give rise to the bulk of the Mgl absorption in the form of very narrow (6 < 1 km s~^) 
components. This phase also gives rise to the observed full Lyman limit break if log Z — —0.8 
for all the clouds. These "Mgl pockets" are quite small, only ~ 100 AU. 

Prom many different techniques, there is conclusive evidence for the existence of small- 
scale structure in the Milky Way ISM down to scales of 100 AU and lower. The techniques 
include high-resolution absorption studies toward globular clusters and binary stars (eg. 
Andrews, Meyer, & Lauroesch (2001), Meyer & Laurocsch (1999), Watson & Meyer (1996), 
and Meyer & Blades (1996)) and Hi 21-cm absorption studies toward extragalactic sources 
and high- velocity pulsars (eg. Faison et al. (1998) and Frail et al. (1994)). The highest- 
resolution absorption studies yielded small b parameters, 1-2 km s^^, for some of the features 
(Andrews, Meyer, & Lauroesch 2001), consistent with our proposed Mgl pockets. 

The observed small, high-density structures have been difficult to reconcile with the 
idea that the ISM is in pressure balance. Their pressures are two orders of magnitude larger 
than the pressure of the warm, ionized inter-cloud medium. In the context of this z = 0.9902 
system, the pressure of these Mgl pockets would be inconsistent with pressure confinement 
by the Mgll clouds. 

There have been several efforts to reconcile the issue of pressure balance in the inter- 
stellar medium. Elmegreen (1997) has investigated a fractal structure model, driven by 
turbulence, that leads to significant clumping on small scales. Heiles (1997) has considered 
modifications of the standard cooling paradigms as well as alternative geometries. Walker 
& Wardle (1998) investigated the idea that the small clouds are self-gravitating. Regardless 
of the mechanism of formation and stability of these small structures, it is clear that they 
are pervasive in the disk of the Milky Way. Therefore, it is not surprising that we would see 
them in absorption through other galaxies such as this one probed by the PG 1634 + 706 line 
of sight. We should note, however, that since there is no definite theory for their existence, 
our assumption that these clouds are in equilibrium may only be an approximation to a more 
complex physical situation. 

A clue about the spatial arrangement of the cold and warm absorbing gas in the z = 
0.9902 system comes from the relative sizes we have derived from the photoionization models. 
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The sizes of the Mgl pockets are 10^-10^ times smaller than those of the Mgll clouds, yet we 
find several of them along this line of sight. This suggests a sheet-like structure such that 
the covering factor is large, or perhaps a fractal structure as proposed by Elmegreen (1997). 
In addition, depending on the geometry, such small Mgl clouds could only partially cover 
the quasar beam, a phenomenon generally associated with absorbing clouds intrinsic to the 
quasar (e.g. Barlow & Sargent (1997), Hamann et al. (1997), Ganguly et al. (1999), and de 
Kool, Korista, & Arav (2002)). The quasar continuum source beam size would be ~ 10-100 
Schwarszchild radii at the distance of the quasar, ~ 20-200 AU for a 10^ Mq black hole. 

In fact, the large 14^(Mgl)/l^(MgIl) ratio that led to our proposal of the existence of the 
Mgl cloud phase is common in strong Mgll systems. In about half the strong Mgll clouds, 
a simultaneous Voigt profile fit to the Mgl and Mgll profiles yielded A'"(Mgl)/A^(Mgll) 
larger than could be reconciled with a single-phase photoionization model (Churchill 1997; 
Churchill, Vogt, & Charlton 2002). Also, Rauch et al. (2002) found a large Mgl column 
density for the strong Mgll system at z = 0.5656 toward Q 2237 + 0305, and found that 
it implied either an extremely large density, or a large iV(Hl). They favored the latter, 
suggesting that this is a strong Lyman limit or damped Lya system. Our results here imply 
that the former solution, of high-density pockets, is also reasonable. 

The existence of these tiny "Mgl pockets" may also provide a hint about the nature 
of damped Lya absorbers, especially those at relatively low redshifts. Unlike strong Mgll 
absorbers, which almost always have a L > 0.05L* galaxy within AOh'^ kpc (Bergeron & 
Boissc 1991; Bergeron et al. 1992; Lc Brun et al. 1993; Stcidel, Dickinson, & Persson 
1994; Steidel 1995; Steidel et al. 1997), DLAs have a variety of types of galaxy hosts (Rao 
& Turnshek 2000), including dwarfs and low- surface- brightness galaxies (LSBGs) (Steidel 
et al. 1997; Turnshek et al. 2001; Kulkarni et al. 2001; Bowen et al. 2001; Bouche 
et al. 2001). If the strongest Hi absorption (the DLAs) occurs in a continuous medium 
in the central region of a galaxy, it seems difficult to explain why Lyman limit absorption 
(17.2 < logA'^(Hl) < 20.3) would not occur in a relatively large annulus surrounding the 
central regions of dwarfs and LSBGs. However, if the DLA is produced in small, high- 
density pockets of material, it is plausible that the surrounding regions could be evacuated 
of gas in galaxies with small potential wells. Thus, dwarfs and LSBGs could give rise to 
DLA absorption but not to a significant fraction of ordinary Lyman limit systems. 

The idea that DLAs come from small pockets is supported by other observations. Lane 
(2000) have resolved narrow components (2-3 km s~^) in an Hi 21-cm absorption profile of 
the DLA at z = 0.0912 toward B 0738+313. With some thermal scaling, it would be expected 
that the associated Mgl components would be even narrower, consistent with what we have 
found in our z — 0.9902 absorber. Furthermore, the DLAs with larger Mgl equivalent widths 
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tend to be the ones with many components in their 21-cm absorption profiles (Rao 2002; 

private communication), and many of these components are also quite narrow (Lane 2000). 
Finally, the high-ionization absorption associated with DLAs is similar to that of classic 
strong Mgll absorbers, as if it is due to parts of the galaxy that are unrelated to the DLA 
region (Churchill et al. 2000). DLAs could simply be small concentrations of gas that exist in 
large numbers in many kinds of environments. Our strong Mgll absorber a.t z = 0.9902 could 
be a less extreme (in its A^(Hl)) and more common version of the same sort of structure. 

5.3. Alternatives to Sepeirate Mgi and Mgii Phcises 

In our favored model, the Mgll phase (with 6 < 6(Mg) < 15 km s^^) severely under- 
produces Mgl, while the Mgl phase (with b ~ 0.75 km s~^) produces only negligible Mgll. 
This might suggest that there would be an alternative model with intermediate b values so 
that Mgl and Mgll can be produced in the same clouds. However, we find that this is not 
the case. Although, with a smaller value of b the observed equivalent width ratio of Mgl to 
Mgll can be matched, the fit to Mgll is unsatisfactory (too deep and too narrow). 

Alternatively, we could relax our assumption that the Mgll profile is fit with the mini- 
mum number of component clouds that produce an adequate fit to the data. A larger number 
of clouds could be used, with smaller b values (perhaps 2-3 km s""^). It would seem that we 
could increase the number of clouds (thus decreasing the maximum A^(Mgll) for a cloud) so 
that Mgll is no longer on the fiat part of its curve of growth. With Mgll on the linear part of 
its curve of growth, M^(Mgl)/H^(Mgll) would be increased. In fact, this is not an acceptable 
solution for this system. As the number of clouds is increased, there is a competing effect 
that reduces W^(Mgl)/W^(MgIl). Increased blending of the superimposed clouds affects the 
Mgl more significantly than the Mgll. This effect dominates so that we are unable to find 
a model with many moderate b clouds to fit the data. 

We conclude that a two-phase model, with separate phases producing the Mgl and Mgll 
absorption, is the simplest suitable model that can fit both the Mgl and Mgll profiles and 
those of other intermediate ionization transitions. 

5.4. Broad Hi phase 

A broad Hi component, with 6 ~ 60 km s~^ was proposed to produce a self-consistent fit 
to the Lya, Ljf3, and Ly7 profiles. In order that this component does not over-produce any 
of the high-ionization metal-line transitions, it must have a low metallicity, log Z < —2. It is 
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consistent with a photoionization model, but the ionization parameter is poorly constrained 
(since this phase does not give rise to metals). The size of this cloud is therefore also poorly 
constrained. If log U — —2.5 the size would be ~ 60 kpc. 

In most other strong Mgll absorbers, the CiV absorption is too strong, and the com- 
ponents too broad, for it to arise in the same phase as the Mgll (Churchill et al. 2000). 
These absorbers also typically require an additional phase to self-consistently fit the Lyman 
series, and it is possible to produce both the high-ionization absorption and the additional 
Hi absorption with the same phase. The need for a broad, high-ionization phase is especially 
evident in cases like the "double" strong Mgll absorber at z = 0.9276 toward PG 1206 + 459, 
in which CiV, Nv, and Ovi absorption are all very strong (Churchill & Charlton 1999; Ding 
et al. 2002). Such strong, broad, high-ionization absorption also characterizes the corona 
that surrounds the disk of the Milky Way (Savage et al. 1997, 2000). 

The broad Hi phase of the z = 0.9902 absorber resembles the Hi component of the 
broad phases proposed for those other strong Mgll absorbers. However, in this case no 
high-ionization metal-line absorption is produced by the broad phase, which led us to the 
conclusion that it has very low metaUicity. With such a low metallicity logZ < —2 as 
compared to its Mgll clouds (logZ ~ —0.8), it seems unlikely that this system is analogous 
to the Milky Way corona. In fact, the low metallicity is more suggestive to a halo structure. 
The b parameter of ~ 60 km s~^ could also be consistent with such an interpretation. 

This leads back to why the system is "C iv-deficient" (and Nv and Ovi-deficient as 
well) . It is more likely that this absorber is produced by a galaxy without a significant corona. 

The alternative, a low-metallicity corona, is less likely since a relatively high metallicity {Z ~ 
—0.8 for the Mgl and Mgll phases) would be indicated for the disk that would produce the 
corona. An early-type galaxy would be consistent with this interpretation, but is certainly 
not a unique solution. It may be only a coincidence that the broad Hi component in this 
CiV-deficient absorber resembles those needed to fit the Lyman series in most other strong 
Mgll absorbers. 

5.5. CoIIisionally Ionized Phase 

Another unusual feature of the z = 0.9902 system is the relatively strong, smooth Siiv 
profile, fit with 6 ~ 17 km s~^. Since the cloud structure apparent in the CiV profile matches 
the Mgll clouds, we infer coUisional ionization with a temperature close to the peak for Si IV 
production, T ~ 60, 000 K. 

It is interesting to note that a similar situation was found in the case of the z = 1.04 
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multi-cloud weak Mgll absorber along this same quasar line of sight (Zonak et al. 2002). In 
that case, the Sim profile was smooth and featureless, and also much stronger than expected 
from photoionization models of the other phases inferred for that system. In that case, a 
somewhat smaller temperature of 40, 000 K provided an adequate fit. 

The temperatures of these additional, coUisionally ionized phases place them on an 
efficient part of the cooling curve, so they would have to be quite common to be detected 
during a relatively brief interval. However, it is reasonable to think that cooling layers 
heated by supcrnovae shocks could provide an appropriate environment. More complex 
models than pure collisional ionization have been considered to explain the ratios of high- 
ionization metal-line transitions in the Milky Way disk (Savage et al. 1997). Certainly, 
more complex processes than pure collisional ionization could be involved in this case as 
well. If strong, smooth profiles for one particular transition are commonly found in strong 
Mgll absorbers, they could be compared in more detail to supernova models and used to 
track the evolution of the ensemble of supernovae remnants in the universe. 
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Fig. 1. — The normalized STIS and HIRES profiles of various key transitions are displayed 
in velocity space. The velocity zero-point is at 2; = 0.9902. The solid curve just above zero 
represents the la error bar. The curve superimposed on the data is a synthetic spectrum for 
one of the best models we have considered (summarized in Table 1). The Mgl, Mgll, and 
Fell profiles were observed with Keck/HIRES, and all the others with HST/STIS. 
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Fig. 2. — Equivalent width is plotted vs. column density, both on logarithmic scales. The 
dashed lines represent curves of growth for Mgll 2796, while the sohd ones represent Mgl 
2853 curves of growth. For both transitions, a set of Doppler parameters, b = 8.0, 4.0, 2.0, 1.0 
and 0.5 km s^^, are indicated. Dotted lines represent typical iV(Mgll) and A^(Mgl) values, 
i.e. those for the cloud at —22 km s~^. From the curve it is shown that for our system Mgl 
is typically on the linear part of the curve of growth, while Mgll is on the flat part. 
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Fig. 3. — Various key metal transitions are plotted in velocity space. The STIS and HIRES 
data are histograms, with model curves superimposed. The dotted curves represent the 
Mgl clouds. The dashed curves show the contribution from the Mgll phase. The dashed- 
dotted curves represent the broad Hi phase. The solid curves show the contribution from 
the coUisional component. 
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Fig. 4. — The Lya, Ly/3, and Ly7 profiles are displayed in velocity space. The solid his- 
tograms represent the STIS data. Other curves show the results for our "best model" (listed 
in Table 1). The model curves represent different phases as indicated in the caption of 
Figure 3. 
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10.2 


Mgna 


1 


0.16 


-3.5 


0.02 


60 


13000 


18.6 


16.8 


10.7 


12.8 


12.1 


11.9 


12.0 


16 


6.5 


6.1 


Mgn4 


12 


0.16 


-3.0 


0.006 


1000 


14000 


19.3 


17.0 


10.6 


13.0 


11.6 


13.1 


13.6 


21 


14.7 


14.5 


Mgns 


24 


0.16 


-2.5 


0.002 


500 


18000 


18.5 


15.7 


9.3 


11.5 


8.9 


12.6 


13.4 


18 


6.3 


5.7 


MgHe 


42 


0.16 


-2.5 


0.002 


500 


18000 


18.5 


15.7 


9.3 


11.5 


8.9 


12.6 


13.4 


17 




3.6 


MgHy 


55 


0.16 


-2.7 


0.003 


50 


16000 


17.7 


15.1 


8.9 


11.0 


9.0 


11.7 


12.4 


17 




3.6 


Broad Hi Phase I 


Hu 





0.003 


-1.5 


0.0002 


60000 


38000 


19.6 


15.5 





7.3 





9.5 


12.8 


58 


52.6 


52.5 


CoLLisiONAL Ionization Phase 


Siivi 


17 


0.16 








63000 


19.2 


15.1 





10.3 


12.0 


13.5 


13.1 


36 


17.2 


16.4 



Note. — 



Column densities are listed in logarithmic units. 



